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Abstract 



An expression is given for the electronic polarization of an insulating crys- 
tal within the ultrasoft pseudopotential scheme. The pseudopotential charge- 
augmentation terms modify the usual Berry-phase expression, and also give 
rise to a second term that takes the form of a conventional expectation value. 
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Typeset using REVTeX 



Several years ago, we presented a formulation that allows for the calculation of the 
electronic polarization of an insulating crystal in the context of a conventional electronic 
band-structure calculation.§'l It was shown that the contribution of a given band to the 
electronic polarization could be expressed in terms of the charge center of its associated 
Wannier function; or equivalently, in terms of the Berry phase of the cell-periodic Bloch 
function as the wavevector k is adiabatically transported by a reciprocal lattice vector. 
Resta has given a review of this theory.! 

The above theory can be applied straightforwardly in the context of an all-electron (e.g, 
augmented plane-wave) scheme, or any norm-conserving pseudopotential scheme. However, 
in the ultrasoft pseudopotential scheme,§^ the pseudo-wavef unctions do not obey a conven- 
tional normalization condition, and the Berry-phase expression must be modified. 

The purpose of this Report is to present the appropriate formulation of the electronic 
polarization within the ultrasoft pseudopotential scheme. The expression has already been 
used in applications,^^ but has never previously been reported. 

Let / = {R, r} label the atomic sites in the crystal, where R is a real-space lattice vector 
and r is the extra fractional translation needed to specify the location of the atom in the 
cell. Note that r will double as a species label for the atom. The basic ingredients of the 
USPP schemei^i are the projectors /3^,i(r) and the charge augmentation functions (5T-ij(r) 
for an atom of species r. Here z is a composite index labeling the angular momentum [Im) 
and radial indices of the projectors. The matrix quantities Qr,ij and d^-^ij are then defined 
via 



Qr,ij = J drQr,ij{r) (1) 

and 

dr,ij = J drrQrAjir) . (2) 

The projectors and charge augmentation functions of atom / located in the crystal at R + r 
are 

/5P(r)=/3.,(r-R-r) (3) 

and 

Q§(v)=Q^,,ir-R-T) . (4) 
The charge density operator in the crystal is thusi 

K(r) = |r)(r| + E5:gS?(r)|A^'))(/5f| (5) 

I ij 

and the number operator is S" = / dr K{r) or 

^ = i + EEQ^..I/5r)d • (6) 



//^(Oi 

I ij 



The Bloch pseudo-wavefunctions ip]^ are written in terms of cell-periodic functions Wk in 
the usual way. 
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7Ak(r) = e^''-'-Mk(r) . (7) 

The normalization {ipk\S\ip\^)y = 1 is relative to the number operator of Eq. (|^). (The 
notation ()^ indicates that the expectation value or matrix element is taken per unit cell of 
volume V.) Substituting Eq. (|]) and using the definition 

/?5)(r) = e-''^-(--)/3.,(r-r) , (8) 
the normalization of the cell-periodic Bloch function becomes 

1 = (wk|«k)v. + EE Qr., («kl/35^) (/35Vk) . (9) 

T ij 

Now, from the Bloch pseudo-wavefunctions, one can construct a localized Wannier 
pseudo-wavef unction in the usual way, 

where the integral is over the Brillouin zone (BZ). It is straightforward to check the normal- 
ization {w\S\w) = 1, viz.. 



{w\s\w) = {w\w) + J2y: Qr., HP?) (/?]'V) 

I ij 



V 



T ij ) 



(2vr) 



|^^rfk|K|nk),. + EEQ^..(«k|/?l?)(/3|?l«k)| . (11) 



The second line results from substituting Eqs. (H) and (0), and the third is obtained from 
the completeness relation 

E e-^C^-'^')-^ = ^ (5(k - k') . (12) 

R ^ 

Using the Bloch normalization of Eq. (|^), it then follows that {w\S\w) = 1. 

The contribution of this band to the electronic polarization is given in terms of the 
Wannier center as {—2e/V){r), assuming spin-paired electrons.Ehi Thus, our remaining task 
is to calculate the location (r) of the Wannier center. Using Eqs. (|1|), (|^), (|), and (|^), 

(r) = J drr {w\K{r)\w) 

= {w\r\w) + J2J2 [iR + r)Qr,j+d,,j] {w\pP) {P?\w) . (13) 

I ij 



3 



Substituting Eq. (p^, one obtains, in analogy with the second hne of Eq. (|TT 



Replacing R+r or R+r by — iVk' acting on the exponential factor, then using an integration 
by parts to cause Vk' to act on the -Uk' and -* instead, and applying the completeness 
relation (0), one finds 



(r) = (r)^^ + (r)^v (15) 



where 



{rf" = EE ("-IfS") (fS'lfk) . (17) 



and 



A(k) = t (MklVkl^k),. +T.T.'Qr,^J («k|/9??) Vk • (18) 



The decomposition (|I^) is into "Berry-phase" (BP) and "expectation-value" (EV) terms. 
The latter involves no coupling between k-points [cf. Eq. (|17D] and is thus obviously gauge- 
invariant, i.e., invariant with respect to a k-dependent change of phase of the Wk. The BP 
term is also gauge- invariant, as can be seen by substituting Wk — > Wk x exp[— i7(k)] into 
(p!8|) and using the Bloch normalization condition (^ to show that A(k) — A(k) + Vk7(k). 
The quantity A(k) thus plays the role of the "Berry connection" or "gauge potential" of the 
Berry-phase theory.i'0 Since exp[— Z7(k)] is periodic in k-space, it follows that 7(k) must 
take the form of R ■ k plus a periodic part (R being a lattice vector), so that (r)^^ is clearly 
invariant modulo a lattice vector. 

In practical calculations, a discrete mesh of k-points is typically used. The expectation- 
value term is easily evaluated on any mesh of A^^ k-points as 

(r)"^ = l^EEE d.., (t^kl/?!?) (/5i?|«k) . (19) 

jVfc T- -J 

Following Ref. |l|, the Berry-phase term is evaluated one component at a time. That is, the 
component of (r)^^ along some primitive reciprocal lattice vector G is evaluated by choosing 
a 2D mesh of A^^^j. strings of k-points running parallel to G, and calculating 

(ril)^^ = ^ E In n {(^kl^k+b). + EE Qr.j («k|4?) (/?5+'Vk+b)| . (20) 

•^^^fc-L k^ I r ij ' } 
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Here ry and k\\ are components of r and k parallel to G, with k = + A;||G, and b being 
the separation between neighboring k-points along the strings. As in Ref. |I], this discrete 
version of the Berry-phase term is manifestly invariant with respect to a change of phase of 
any one of the u^, since each appears once in a bra and once in a ket. 

For the case of multiple discrete bands, the final expression for the polarization is 

P = Pion-^E[(r)f + (r)r] , (21) 

where the sum is over occupied bands n. For the case of compound multiple bands that 
have degeneracies at certain locations within the Brillouin zone, the EV terms present no 
difficulties, but the BP evaluation requires care. One can again follow the example of Ref. |l| 
to express the Berry-phase contribution to the electronic polarization as 

P BP = ^ E Im In n det M^^) (22) 

where 

= («™,k|Un,k+b).^ + E E Qr,, («™k|/5S^) (4';+^Vn,k+b) (23) 

T ij 

and the indices m and n run over occupied bands. 

In summary, the formulas for electronic polarization given in Ref. have been general- 
ized to the case in which ultrasoft pseudopotentials have been used in the band-structure 
calculation. The resulting formulas are easily implemented in practice. 
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